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Abstract 

Reconstructing the evolution history of the dark energy equation of state pa¬ 
rameter w{z) directly from observational data is highly valuable in cosmology, 
since it contains substantial clues in understanding the nature of the acceler¬ 
ated expansion of the Universe. Many works have focused on reconstructing 
w{z) using Type la supernova data, however, only a few studies pay attention 
to Hubble parameter data. In the present work, we explore the merit of Hubble 
parameter data and make an attempt to reconstruct w{z) from them through 
the principle component analysis approach. We find that current Hubble pa¬ 
rameter data perform well in reconstructing wlz); though, when compared to 
supernova data, the data are scant and their quality is worse. Both ACDM and 
evolving w{z) models can be constrained within 10% at redshifts z < 1.5 and 
even 5% at redshifts 0.1 < z < 1 by using simulated H{z) data of observational 
quality. 
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1. Introduction 


There have been a handful of works reconstructing the dark energy equation 
of state parameter w{z), defined as w{z) = P/p, where P and p represent the 
pressure and energy density, from observed distance modulus data (e.g., Type 
la supernova (SN la) data). Basically, these works can be divided into two cat¬ 
egories: (i) Assume a particular model about w{z) so that the model quantities 
(e.g., luminosity distance cLl) can be compared directly with the observations, 
(ii) Characterize the data with an underlying model and then reconstruct w{z) 
numerically. For category (i), the assumed model can be expressed as a pa¬ 
rameterized form of a basis representation or a function distribution represen¬ 
tation. The parametric forms are currently the most commonly used, which 
either assume w{z) is a constant or allow for a redshift dependence, such as 


(z) = wq — wiz/{\ z), where wq and w\ are constants 


QEi- 


Hojjati et 


al. [3| used wavelets to study the redshift evolution of w{z). In that work the 
Principal Component Analysis (PCA) technique, which has the advantage of 
detecting local features such as bumps in w(z), was applied to decorrelate the 
wavelet coefficients. Huterer et al. [ 4 I presented a piecewise constant description 
of w{z) by binning redshifts, and estimated the corresponding band power. The 
PCA was also used to make the band power uncorrelated so that the eigenmodes 


■ B 


can be written as linear combinations of redshift bins. Holsclaw et al 
troduced a non-parametric reconstruction method by modeling w{z) through a 
Gaussian process (GP) and estimating the GP hyperparameters with a Bayesian 


inference and Markov chain Monte Carlo sampling. Zhao et al. [gl combined 
GP and redshift binning method to feature w{z), in which the w value in each 
bin is treated as free parameter to be determined. To make w{z) differentiable 
for the calculation of the dark energy perturbations, they interpolated between 
the bins with narrow tanh functions.For category (ii), Daly et al. B devel¬ 
oped a simple numerical method to determine w{z) from the data, which made 
no assumptions about the underlying cosmological models. Their approach is 
model-independent. However, the final results are noisier and highly sensitive to 
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the amount and quality of the available data. To suppress the noise, Shafieloo 
et al. Q smoothed the SN la data by applying a Gaussian kernel over the whole 
redshift range before performing reconstruction. Majumdar et al. Q expressed 
the distance modulus in a parametric form, fi = Slog [z{l + az)/{l + bz)] + 
and used the standard minimization to determine the parameters (a, 6, ^q). 
Clarkson et al. [l^ fitted the data before reconstruction using eigenmodes which 
were transformed from a set of primal basis functions with PC A. 

With observational data getting richer and more diverse, recently researchers 
have begun to study dark energy reconstruction with multiple sets of measure- 


n 


to 


ments. Holsclaw et al. nil extended their Gaussian process method 
include a diverse set of measurements: baryon acoustic oscillations (BAO), cos¬ 
mic microwave background measurements (CMB), and supernova data (SNe), 
while observational Hubble parameter H{z) data (OHD) were used by other 
researchers Q, [d , 121. 


So far, most research has focused on observational data acquired from SNe, 
CMB, and BAO, especially from SNe, or their hypothetical counterparts. On 
the contrary, the problem of reconstructing dark energy from Hubble parameter 
data has attracted less attention of the cosmology community. How valuable are 
OHD in modeling the Universe? In fact, some pioneers have made explorative 


works in this area. Ma et al 


• 0 


studied the viability of constraining the 


cosmological parameters using OHD and concluded that as many as 64 data 


■ 0 


points are sufficient to get a result comparative to SNe data. Pan et al. 
discovered that to some extent the OHD were better poised to give information 
about the nature of dark energy than the SN la data. Encouraged by these 
pioneering works, in this paper we attempt to reconstruct dark energy solely 
with the help of OHD. 

The paper is structured as follows. In section[2J we analyze the advantage of 
OHD. The reconstruction method we adopt is outlined in section [3] Our results 
are presented in sectional and finally we make some conclusions in section [5] 
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2. Merit of the Hubble parameter 


Although, as previously mentioned, some authors have considered the role 
of OHD in reconstructing dark energy with combined data sets, few researchers 
have tried to use OHD solely to undertake this task [1^. This is partially 


due to the smaller number and lower quality of OHD compared to the dis¬ 
tance module data. In contrast to the 580 SN la data released by the su¬ 


pernova cosmology project 


H, 


to date, the most comprehensive sample of 


OHD contains only 37 measurements (see Figs. [T|). They are collected from 
[l2) 18, 3, 3, 21,1^, 3, Q, H, Q) most of which are compiled in 261. Cur¬ 
rent OHD are obtained primarily by the cosmic chronometers method (30 mea¬ 


surements) [D 
to extract H{z) [2l|, 


15j. Radial BAO peaks detection is another method 
23j. The cosmic chronometers method is affected by 
the systematic biases in the age determination of galaxies, while for the radial 
BAO method, the systematics stem from various distortion effects subjected 
to the autocorrelation function. Since a fiducial flat Lambda cold dark matter 
(ACDM) model and its parameters are used to convert redshifts into distances 


and to gauge the comoving BAO scales in the selected redshift slice 


27|, the 


radial BAO method is essentially model-dependent, and the possible circular 
logic issue also contributes to the systematic errors in this method 2^. Since 


OHD measurements coming from different techniques may suffer from different 
systematic effects, the 30 measurements that come from cosmic chronometers 
are adopted in our analysis (see FigH]). They all have been compiled in Ref. 

H- 

OHD (with mean relative error of 19.8%) have generally larger error bars 
than SN la data (with mean relative error of merely 0.55%). However, OHD 
H{z) have more direct and efficient power to constrain w(z) compared to SN la 
dL{z) data. w{z) is written in terms of dL{z) as Equation ((T|) |10l |: 


uj(z) = {2{1 + z){l + nkD^)D"-[{1 +z)^nkD'^ 

+2(1 + z)nuDD' - 3(1 + DfcD2)]D'}/{3{(l + zf 
X [flk + (1 + z)nm]D'^ - (1 + nkD^)}D'}, (1) 
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Figure 1: OHD data and the best-fit ACDM model. The black solid curve represents 

the ACDM model best-fitting the data, which gives Qm = 0.36, Qk = —0.13 and Hq = 
69.1 km/s/Mpc. 

Here ^Irn flk normalized density parameters of matter and curvature 
and prime denotes the derivative with respect to z. D(z) = dL(z)Ho/(c(l + z)) 
is the normalized comoving distance and Hq = H{0) is the Hubble constant. 
Note that the second derivative of D{z) is a prerequisite in Eq. ([T]) even in 
the case of Hfc = 0. We can decompose Eq. © into two simpler equations by 
introducing the Hubble parameter H{z) as the intermediate quantity, which can 
be expressed as Equation (l2|), 

H{z) = (2) 

and Equation Q, 

12HH'{l + z)-3H^ + H§nkil + z)^ 

Eq. ([3]) shows that to calculate w{z) from H{z), only the hrst derivative of H{z) 
is needed. The derivative operator functions as a high-pass filter and can thus 
heavily amplify the observational errors. Therefore, to reconstruct w(z), H{z) is 
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^referable to (Il^z) when regarding the suppression of the propagation of noise 


29|. Such an advantage of H{z) makes it a promising method to reconstruct 


w{z) directly from H{z) instead of dL(z), although the current quality of H{z 


data is comparatively poor. It should be noted that Pan et al. 


14| have given 


an estimate about the advantage of the first derivative of data over the sec¬ 
ond derivative in reconstructing dark energy. They studied the efficacy of four 
different parameters: the deceleration parameter, the equation of state of dark 
energy, the dark energy density, and the geometrical parameters Om{z), which 
can be expressed as Equation O 0: 




( 4 ) 


In discriminating theoretical models of dark energy using supernova data, they 
found that the geometrical parameters Om{z), a cosmological parameter con¬ 
structed from the first derivative of the data, for which the theoretical models 
of dark energy are sufficiently distant from each other, performs best in recon¬ 
structing dark energy from SNe data. 


3. Methodology 

Most existing reconstruction approaches making use of the observed distance 
modulus data can also be applied to OHD. Here we apply the framework of 
weighted least square (WLS) and PC A Q. Given a set of H{z) data X = 
{xi,X 2 , , we fit them with a smooth analytical function, and then use 

it and its first derivative to reconstruct w{z) according to Eq. Assume the 
covariance matrix C of the data X satisfies Ca = af for i = 1,2,..., AT and 
Cij = 0 for i ^ j. We start with a set of N primary basis functions Pn{z), 
and their summation o,nPn{z) is used to fit X through the weighted least 

square technique, with the weight matrix being C~^. Minimizing the weighted 
squared residual TZ = {X — PA)'^C~^{X — PA) by solving dTZjdA = 0, we get 
the coefficient vector A = (oi, 02 ,..., ajv)^ as shown in Equation ((S]): 

A = (P'^C-^P)-^P'^C-^X, (5) 
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where the K x N measurement matrix P = {pi,p 2 , ■■■,pn) is composed of the 
basis function set. Meanwhile, according to the least squares theory, the inverse 
covariance matrix of coefficient vector A is given by Equation ([51): 

Cov-'^{A) = P'^C-^P. ( 6 ) 


Then the PCA is implemented by calculating eigenvalues and eigenvectors 
of Cov~^{A): Cov~^{A) = EAE'^, where E = (ei, 62 ,..., bat) is the matrix 
of eigenvectors and A the diagonal matrix with diagonal entities representing 
eigenvalues with decreasing order. A set of N optimal basis functions Q = 
( 91 , 92 , 9 Ar) are produced by transforming the primary functions through E: 
Q = PE. That is given by Equation ([7]), 

N 

Qnjz) = emnPm{z),n = (7) 

m—1 

where Cmn is the mth component of the nth eigenvector. The first M {M < N) 
eigenfuncitons Q' = ( 91 , 92 , 9m) are usually of interest since they are the 

principal components and represent the main features of the data. We use them 
to ht X once more with WLS and get new coefficients A' = (a^, 02 ,..., 
given by Equation ([8]) : 


A' = {Q'^C-^Q')-^Q'^C-^X. 


( 8 ) 


Then the function ^nQri{z) gives the analytic representative of the data 

X. Thus, each choice {N,M} will give a particular reconstruction of w{z). 

To balance the risk of getting a wr ong 'w{z) and the risk of over-fitting. 


we adopted the formalism discussed in lid- That is, a combined information 
criterion is used to cull out the most suitable combinations of {A^, M} with 
N G [2,10] and M G [2,iV], which is defined as Equation ([S]): 


CIC = (1 - s)AIC -f sBIC, 


(9) 


where AIC is the Akaike information criteria and BIC the Bayesian criteria. 
This framework does not single out one particular combination of {N, M} that 
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possesses minimum CIC, but selects a family of M} reconstructions which 
satisfy Equation 

CIC < ClCniin + K, (10) 


where k is a constant being set to 5 through the paper. From the selected family 
of {N,M} reconstructions the averaged w{z) and 1 —cr errors are inferred. Thus, 
K controls the size of the family of {N,M} combinations and s gives the weight 
for the importance of the Bayesian criteria relative to the Akaike criteria. 

Both OHD (denote by Ho{z) hereafter) and simulated H{z) data (denote 
Hs{z) hereafter) are tested in our experiments. Moresco 2^ performed a H{z) 
simulation of a Euclid-like survey to forecast the expected improvement on es¬ 
timating cosmological parameters with future data. Here, we get Hs{z) based 
on currently observed data by using the method proposed in [l2|. We assume 
Q,ra = 0.30 and = 0 and the Hubble constant Hq is set to the newest 
value of 67.4km/s/Mpc 25|. In this method, we take the relative error <jh/H 
as a random variable which satisfies the z-independent Nakagami distribution 
/m (x;m, H). By htting OHD, we get the following distribution parameters: 
m = 0.63 and H = 0.038. The final Hs{z) data should follow the Gaussian dis¬ 
tribution Hs{z) ^ N{H^{z), cfh): where is the fiducial H which gives the 

simulated data. In total, 400 data points with 2 evenly distributed in the range 
of [0, 2.0] are generated. The uncertainties of the Hs{z) increase with z. As a 
result, the simulated data set is qualitatively comparable with observations. 


4. Results 

4.I. Results on simulated data 

We hrst test the power of the Hs{z) data with the above technique. The 
primary basis functions are Pn{z) = [■^/(l + 2 )]"“^, n = I ^ A^. Other types 
of functions, such as and 1/(1 -I- 2 )"“^, are also tested but led to similar 
results. The function set 2 ;”“^ leads to slightly larger reconstruction errors, 
while the reconstruction errors given by [1/(1 + z)]’^~^ are slightly smaller than 
[z/{l + z)Y^~^. In FiglJl we show the first 10 normalized eigenfunctions obtained 






Figure 2: The first 10 eigenfunctions for the Hs{z) data, with N = 10. Note the numbers 
aligning the ordinate mark the order of these eigenfunctions other than their own values. 

by fitting the Hs{z) data with PCA, with the underlying model being w = -1 
and N = 10. We find that, for a fixed N, the shapes of eigenfunctions vary from 
smooth to oscillatory. The fiducial model used in this procedure has the same 
parameters as those used in the simulation procedures, i.e., Hq} = 

{0.30,0.0,67.4} and Wflduciai = —1- We set /t = 5 and study three cases s = 
0.1, 0.5, and 1.0, respectively. FigI3] shows the reconstructed w{z) (denote by 
whs{z) hereafter). On the whole, the results are quite good for all s values. 
Fig l3] shows that s takes us from conservative models where s = 1 to more wild 
models where s = 0.1. This is due to the fact that models minimizing the BIC 
criteria tend to become smooth with tight error bars, while models with low 
AlC values are oscillatory and have large errors. 

It is expected that the number of data points should have a crucial impact 
on the quality of the reconstruction. As shown in FigUl the reconstruction 
error drops drastically with an increasing number of the simulated data points, 
except for the number greater than 200 or so. For 200 or more data points, the 
resulting errors are comparable to what are shown in FigjSl We can see that the 
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Figure 3: Reconstruction of w{z) using Hs(z) data. The solid curves represent the recon¬ 
structed w{z) and the dashed curves represent 1-cr errors. The gray horizontal line denotes 
w = —1 

currently observed 30 Ho{z) data points, with current measure quality, are really 
too sparse. We expect more OHD data points to be measured in future. On the 
other hand, 400 data points seem to be enough for improving the reconstruction 
quality, and further increasing the number of data points will merely lead to a 
marginal improvement. Further improvement on the reconstruction result needs 
further improvement on data quality. 

However, it is unlikely that future surveys will provide hundreds of H{z) 
measurements. On the contrary, current measurements are basically systematic- 
dominated, so there is space for improvement in the direction of getting smaller 
errorbars, even at higher-z. Thus, we give the H{z) simulation by keeping the 
number 30 of data points fixed but reducing Hs{z) errors. How the recon¬ 
structed w(z) errors vary with different Hs(z) errors is shown in Figlb] We can 
see that, to make the reconstructed w(z) error less than 0.05 at 2 ; = 1, which is 
equivalent to the reconstruction accuracy obtained by 400 Hs(z) points of cur¬ 
rent measurement quality (see FigHJ, the relative error of Hs(z) should achieve 
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Figure 4: The reconstruction errors of w(z) at z = 1. Here w(z) are reconstructed by Hs{z) 
data with different number of data points. We set s = 1.0 and k = 5. 


4%. Obviously this is a big challenge for future surveys. We alternatively test 
50 Hs{z) points and find the needed relative Hs{z) error should be 6%, and for 
100 Hs{z) points the relative error 12% is necessary. Considering the feasibility 
of future implement, 100 OHD measurements with quality of 12% are desirable. 

For comparison a hypothetical data set of distance modulus fi{z) data is 
also constructed [l^, which consists of 2000 data points evenly distributed in 
the redshift range 0.08 < z < 1.7, and 300 data points in range [0.03, 0.08] 
The uncertainties are composed of a constant statistical error Umag = 0.15 
and systematic errors linearly drifting from (Tsys = 0 — 0.02. The primary basis 
functions Pn{z) = [z/{l + z)Y‘~^, n = 1 N, are transformed again to a set of 
eigenmodes with PCA, which are suitably adapted to the simulated /r(z) data. 
The cosmological parameters Hq} = {0.3,0.0,65.0} are assumed, ac¬ 


cording to 10|. Different assumptions of Hq in H{z) and ix{z) simulations (Tfg is 
set to 67.4 and 65 respectively) don’t affect the results. Fig [5] demonstrates the 
comparison between w{z) errors reconstructed from Hs{z) data and those recon¬ 
structed from simulated p{z) data. Red curves represent i7-reconstructed errors 
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Figure 5: The reconstruction errors of w(z) at z = 1. Here w(z) are reconstructed by Hs{z) 
data with different quality but fixed 30 data points. We set s = 1.0 and k = 5. 

and blue curves represent //-reconstructed errors. We find that for moderate 2 
values {z ^ 0.4 — 0.8), both i?-reconstructed errors and ^-reconstructed errors 
are relatively low. When z is greater than 0.6, however, the i/-reconstructed 
errors are lower than /^-reconstructed ones, and with the increase of z, the dif¬ 
ferences between them get more remarkable. Increasing s also vastly magnifies 
the differences. For smaller s, in contrast to the iJ-reconstructed errors, the 
/^-reconstructed errors will grow uncontrollably as z increases. Only in the low- 
z region (z < 0.4) does Hs{z) data produce whs{z) errors greater than /i(z) 
data. In addition, both /r-reconstructed and iJ-reconstructed errors grow with 
decreased z in the low-z region (z < 0.4), and the differences between them are 
acceptable. Taken together, Figl6] (see also FigH]) shows that we can constrain 
the ACDM model within 10% at redshifts z < 1.5 and even 5% at O.I < z < 
I by using the simulated H{z) data. Thus, we conclude that using Hs{z) data 
can get better results than that using simulated fj,(z) data. Note that the mean 
relative errors of Hs{z) data are about 0.14 ~ 0.16, while the relative errors of 
/r(z) data are about 0.003 ~ 0.004, indicating the former is nearly two orders 
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Figure 6: Errors of w(z) reconstructed from simulated data. Red curves correspond to Hs{z) 
data and blue curves correspond to fJ,{z) data. The solid, dashed, and dotted lines denote s 
= 1.0, 0.5, and 0.1 respectively. 


larger than the latter (we compare the relative errors rather than the absolute 
values since Hs{z) and fi{z) are different physical quantities). Hence, the com¬ 
parison convincingly validates the merit of Hubble parameter data, as stated in 
section [2l in reconstructing the dark energy equation of state. 

Based on the analysis in section [2l the reconstruction of w(z) from /i(z) 
data can be split into two procedures: hrst reconstructing II(z) (denote by 
Hr{z) hereafter) from fi{z) using Eq.Q, and then deriving w{z) from Hr{z) 
using Eq.(l3]). One may imagine that the error propagation and amplihcation 
taken in the first procedure will induce errors of Hr{z) similar in value to, or at 
least at the same order as, errors of Hs{z) data, so they would produce com¬ 
parable errors (see Eigl6]) for the resulting w{z). However, EiglT] shows this is 
not the case. We find from this figure that, compared to the errors of Hs{z) 
data, Hr{z) errors are quite small. We give the explanation as follows: hrst note 
that both Hs{z) and Hr{z) data are noisy, with errors measuring the extent of 
noise (see FigiT]). Because we use the hrst M < N eigenfunctions, which en¬ 
compass the dominant features in the data, to ht the data, and throw away the 
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redshift z 


Figure 7: Errors of Hs{z) data (red) and errors of Hr{z) data (blue). The solid, dashed, and 
dotted lines denote s = 0.1, 0.5, and 1.0 respectively. Note that the errors produced by H{z) 
simulation are random numbers. 


higher ones which contain noise-induced oscillations (see Methodology), noises 
are suppressed in the resultant whs{z). However, this is the case for whs{z) 
reconstruction from Hs{z) data but not the case for w^{z) reconstruction from 
Hr{z) data. By resampling the parameters in ^{z) data fitting, noises are to a 
large extent suppressed in the first procedure of w{z) reconstruction from /i(z) 
data (i.e., the Ijl{z) — Hr(z) procedure), resulting in relatively low Hr{z) er¬ 
rors (see FiglT]), while the second procedure (i.e., the Hr{z) w{z) procedure) 
essentially maintains all the noise-induced features inherited from the first pro¬ 
cedure. In a word, it is the noise-suppression effect in the reconstruction of 
w{z) from Hs{z) data, which is lacking in the second Hr{z) w{z) procedure 
of w{z) reconstruction from ^(z) data, that makes the w{z) reconstructed from 
Hs{z) behave better, even though the quality of Hs{z) data is much worse than 
Hr{z) data. 

Then we use Hs{z) data to reconstruct two types of evolving Weiz). One 
is a standard slow evolution, expressed as w; = ^[—1 + er/(ln^^ ~1)]; the other 
evolves more drastically, modeled by w = —1 -I- 0.31 sin[12 ln^^“''^^]. The recon- 
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struction result is shown in Figs. [5]and[ni We set s = 0.2 and = 

{0.30,0.0,67.4}. We again see that Hs{z) data outperform, in most cases, the 
simulated SN la data. Similar to the result obtained from the ACDM model, 
the two evolving We{z) models are also constrained within 10% for most red- 
shifts in the range z < 1.5. For the second type of We(z), not only the errors 
constructed by the former are much lower than the latter, but also the former 
constructed We(z) is more faithful to the underlying We{z). It seems that Hs{z) 
data perform better for more oscillatory evolving dark energy. 


^.2. Results on observational data 

Finally we test the performance of OHD. See section 2 for the description 
of the OHD dataset of 30 data points. To make a comparison, the data set 
of distance modulus /i(z) is also acquired. The last distance modulus mea¬ 
surements, which consist of 580 observed points, were produced by Supernova 


Cosmology ProjecJ}! [l^. The parameters of the best-fit fiducial ACDM model 
for OHD data are Hq} = {0.36,-0.13,69.1}, while the parameters 

for htting fj.{z) data are Ho} = {0.32,—0.09,65. 3}[l^. Interestingly, 

the most likely w{z) reconstruction, denoted who{z), favors a transition from 
who{z) < — I at low redshift to who{z) > — I at higher redshift, a behavior 
that is consistent with the quintom model which allows who{z) to cross -1 321. 


The comparison of I — ct errors between iJ-reconstructed and /^-reconstructed 
w{z) is given in Fig ITT] We can see from Fig Hi] that the quality of two recon¬ 
structions are comparable through most of the z region, though in low redshift 
region iT-reconstructed errors are noticeably larger than //-reconstructed errors. 
Roughly speaking, the reconstruction from Ho{z) performs better than the re¬ 
construction from //(z) for relatively large z values. Concerning the much worse 
quality and much fewer observational data points of Ho{z) compared to //(z), 
this result is impressive. 

It is shown from FigHOlthat the shapes of the curves s = 0.1 and s = 0.5 


^ http;//supernova.Ibl.gov/Union/ 
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Figure 8: Reconstruction of evolving dark energy of the first type using Hs{z) data (a) and 
error comparison (b). In (b), the red line denotes the errors constructed by Hs(z) data and 
the blue line corresponds to simulated SN la data. 
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Figure 9: Reconstruction of evolving dark energy of the second type using Hs{z) data (a) and 
error comparison (b). In (b), the red line denotes the errors constructed by Hs{z) data and 
the blue line corresponds to simulated SN la data. 
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Figure 10; Reconstruction of w using Ho(z) data. The solid curves represent reconstructed 
w(z) and the dashed curves represent l-cr errors. 

are quite similar. Accordingly, as shown in Fig llll the difference of error curves 
of who{z) for s = 0.1 and s = 0.5 are small, and the s = 1.0 curve differs 
from the other two iJ-reconstructed curves at a lower extent compared to the 
/i-reconstruction counterparts, especially in high redshift region. This behavior 
reveals that the reconstructions obtained from Ho{z) data are more “stable” 
than those from ^{z) in the sense that they are less dependent on s values. 
Here we can give an intuitive interpretation for it. The parameter s in fact 
controls M values and the number of {N^M} combinations. The smaller s 
is, the higher the probability that the Ms larger, and the bigger the set of 
{N, M} combinations is. Therefore, smaller s results in the emergence of larger 
M in the set of {A^, M} combinations, which leads to more precise data fitting 
and more noisy features reserved. Since ^ is more ’’distant” from w than H, 
even small uncertainties included in the SN la data htting curves may induce 
comparatively large oscillations of the resultant 'Wfj,{z). On the contrary, for 
the W(^z) reconstruction from OHD, the noisy features entering the eigenmodes 
due to the dropped s will not cause the resultant w{z) to oscillate signihcantly. 
Hence, from this aspect the merit of OHD in reconstructing dark energy equation 
of state is established more solidly. 
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Figure 11: Errors of w{z) reconstructed from observational data. Red curves correspond to H 
data and blue curves correspond to /i data. The solid, dashed, and dotted lines denote s = 
1.0, 0.5, and 0.1 respectively. 

5. Conclusions 

We made an attempt in this paper to reconstruct dark energy equation 
of state parameter w{z) from Hubble parameter data within the framework 
proposed by Clarkson & Zuncke Q . We first made a comparison between the 
simulated Hubble parameter Hs{z) data and distance module fj,{z) data. The 
results showed that Hs{z) data can produce more conservative w{z) than the 
simulated ^{z) data. Then the performance of OHD was evaluated against SNe 
data. We found that the errors of the reconstructed w{z), using OHD and /i(z) 
data, are comparable. We also concluded that the number of data points could 
impact the quality of the reconstruction. It turns out that 400 data points 
are sufficient to improve the reconstruction quality using OHD only. Further 
improvement requires more accurate data measurements. 

We have shown that the Hubble parameter dataset alone is potentially ca¬ 
pable of reconstructing the dark energy equation of state. As discussed above, 
to obtain satisfactory w{z) reconstruction, efforts on both the quality and the 


19 






number of OHD are requisite. There are several ways to enhance the quality 
of OHD, of which deeper-redshift, more-complete-sky-coverage LRG survey and 
spectroscopic observations of those identified LRGs give remarkable improve¬ 
ment. It is expected that future surveys will be able to offer a large sample of 
massive and passive galaxies to be used in age-dating. The precise determination 
of the absolute galaxy ages is a quite complicated process, however, that issue 
can be ignored for H{z) measurements because one just needs to take the differ¬ 
ence between the ages in narrow redshift bins so that the systematic bias in the 
absolute ages cancel out. The ongoing GMB observation program, the Atacama 
Cosmology Telescope (ACT) [33|, will produce a catalogue of galaxy cluster can¬ 


didates to be validated up to z ~ 1.5 via the Sunyaev-Zel’dovich effect 


M 


This ensures a promising dataset outperforming current SNIa datasets. For the 
range 1.5 < z < 2, forthcoming surveys such as Euclid [36| and WFIRST [37|, 


will greatly improve the present-day statistics, increasing the number of passive 
galaxies at these redshifts by at least an order of magnitude!^, which will re¬ 
sult in observing approximatively one thousand spectra of passive galaxies at 


these redshifts 


241 . Combining these future observations, it is reasonable to ex¬ 


pect that the Hubble parameter data will be helpful in exploring the expansion 
history of the Universe. 

It should be noted that the results and analyses in this paper are method- 
specific, since we consider just one reconstruction method. Whether similar 
results and conclusions can be obtained by using other reconstruction methods. 


such as Gaussian process modeling Q, is still unknown. 
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